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Abstract. - We use a Ginzburg-Landau free energy functional to investigate diblock copoly- 
mer morphologies when the copolymer melt interacts with one surface or is confined between 
two chemically patterned surfaces. For temperatures above the order-disorder transition a 
complete linear response description of the copolymer melt is given, in terms of an arbitrary 
two-dimensional surface pattern. The appearance of order in the direction parallel to the sur- 
face is found as a result of the order in the perpendicular direction. Below the order-disorder 
transition and in a thin-film geometry, our procedure enables an analytic calculation of dis- 
torted perpendicular and tilted lamellar phases in the presence of uniform or striped surface 
fields. 



Block copolymers (BCP) are macromolecules made up of two or more chemically distinct 
subunits, or blocks, covalently bonded together. The usual (macro-) phase separation occur- 
ring for two immiscible species is not possible for BCP because of the covalent bond between 
the blocks. The most studied BCP are the diblocks, for which the phase diagram is quite 
well understood and was found to consist of disordered, lamellar, hexagonal and cubic 

micro-phases with characteristic length scale usually in the range of dozens of nanometers. 
The prevailing morphologies depend on three parameters: the Flory parameter character- 
izing the incompatibility between the two blocks, the polymerization degree of the chain N, 
and the fraction / = Na/N of the A block in a chain of A^ = Na + Nb monomers. For 
high enough temperatures the system is found in its disordered state. For symmetric melts 
(/ = 1/2) lowering the temperature (or equivalently, raising x) results in a phase transition 
to the lamellar phase through the order-disorder transition (ODT) point. For asymmetric 
melts (/ ^ 1/2) the ordered phases can have, in addition, hexagonal (cylindrical) or cubic 
symmetries 

The technological importance of copolymer thin-film is prominent in diverse fields, such as 
in fabrication of nanolithographic templates [Q, waveguides, anti-reflection coating for optical 
surfaces |^ and dielectric mirrors Q. Therefore, it is highly desired to acquire a better 
understanding of the copolymer behavior in the presence of chemically patterned surfaces. 
Theoretical and experimental ||2l|-p5[| studies have been carried out in order to explore 

how the BCP film thickness, the range and strength of the uniform surface interactions, as well 
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as the bulk parameters x, N and /, compete to produce a rich interfacial behavior. Below the 
ODT, these studies [|l^,^,|l9| revealed parallel, perpendicular and mixed lamellar phases, the 
apparent morphology being the one that optimizes interaction with the surface and periodicity 
frustration. Alignment of the lamellae can also be achieved through the coupling of electric 
fields to local dielectric constant variations , or by shear |27| . 

Let us consider first a BCP melt in its disordered phase (above the bulk ODT temperature) 
and confined by one or two flat, chemically patterned surfaces. Although the bulk BCP is 
disordered above the ODT, there is an oscillatory decay of the A-B correlations and surface 
induced ordering is quite complex. In the vicinity of the ODT this ordering can become 
long range leading to a strong effect. Defining the order parameter (/)(r) = (pAi'i') — / as the 
deviation of the A monomer local concentration from its average, the free energy (in units of 
the thermal energy fcsT) can be written as: 

do = 27r/(7o is the fundamental periodicity in the system, and is expressed by the polymer 
radius of gyration Rg, through qq ~ 1.95/Rg. In addition, r = 2pN {xc — x): — l-5pc'^-Rg/'7o 
and /i is the chemical potential. The second length scale in the system is determined by the 
ratio of two parameters, (r//i)~^^^ ~ {Nxc — A^x)~^^^j and it characterizes the decay of 
surface induced modulations. The Flory parameter x measures the distance from the ODT 
point, having the value Xc — 10.49/iV. Finally, p ~ is the chain density per unit 

volume, and c and u/p are dimensionless constants of order unity Q. 

This and similar types of free energy have been used successfully to describe spatially 
modulated phases p8[-p0|, with applications to amphiphilic systems [^ll^, diblock copoly- 
mers Langmuir films and magnetic (garnet) films [|36| . The free energy, 
Eq. (np^ describes a system in the disordered phase having a uniform (j) — for x < Xc; while 
for X > Xc (but in the vicinity of Xc), the system is in the lamellar phase and is described 
approximately by a single q-niode (j) = (j)q exp{iqQ ■ r). The use of Eq. (|) limits us to a region 
of the phase diagram close enough to the critical point where the expansion in powers of 
(j) and its derivatives is valid, but not too close to it, because then critical fluctuations are 
important 

The presence of chemically heterogeneous surfaces is modeled by adding a short-range 
surface interactions to the free energy, 

F,= J dV, (c7(r,)0(r,) + r,02(r,)) (2) 

where the vector r = defines the position of the confining surfaces. The surface field 
cr(rs) has an arbitrary but fixed spatial variation and is coupled linearly to the BCP surface 
concentration 4>{rs)- Preferential adsorption of the A block is modeled by a constant a < 
surface field, resulting in parallel-oriented layers (a perpendicular orientation of the chains). 
Control over the magnitude of this surface field can be achieved by coating the substrate with 
carefully prepared random copolymers p^,p5t. If the pattern is spatially modulated, cr{rs) ^ 
0, then the A and B blocks are attracted to different regions of the surface. The coefficient 
of the 0^ term in Eq. (^ is taken to be a constant surface correction to the Flory parameter 
X [0 • A positive coefficient corresponds to a suppression of surface segregation of the 
A and B monomers. 

We consider first the semi-infinite system of polymer melt in contact with a single flat 
surface given by = {x,y = 0, z). The ordering effect is expressed in terms of the correction 
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to the order parameter 5(j){r) — (j) — about its bulk value 0^. This correction vanishes in 
the bulk, liniy^oo Scf) = 0, because of the finite range of surface interactions. Moreover, the 
correction should preserve the condition of fixed A/B ratio, namely / 5(j){r)d'^r = 0. The free 
energy can be expanded about the bulk value F[(j)b] to second order in 6(p: F = F[(j)b] + AF, 
with 



AF 



1. 



2 '-V2)J0]'|d3r 
{cr(x, z)5(j)s + Ts {2(l)bS(f)s + 5(t>^)}dxdz 



(3) 

can 

qz), 



54>s = S4i{x, 0, z) is the surface value of 6<p. The surface chemical pattern a{rs) — (t[x, z) 
be decomposed in terms of its g-modes a(x, z) — J2q exp[i [qxX + qzz)], where q = (q^, 
and (Tq is the mode amplitude. Similarly, the correction field is S(f>{x, y, z) — 
^^5(l)q{y) eyipli {qxX + qzz)], and is substituted in Eq. (||). We will first consider a system 
found above the ODT temperature, in the disordered phase. In this case ipb — 0, and the 
X and z integration of the free energy, Eq. (j^), is carried out explicitly. Then, applying a 
variational principle with respect to 54>q results in the following differential equation: 







(4) 



This ordinary differential equation is linear and of fourth order. In the semi-infinite geometry, 
y > 0, the solution to Eq. (^) has an exponential form (50q(j/) = Aq exp{—kqy)+Bq exp(— 
where kq is given by 



fcq — 



q 




Nx) 



1/2 



(5) 



Out of the four complex roots for kq, the two with Re(fcg) < diverge at y ^ cx) and are 
discarded. From the requirement of fixed A/B ratio it follows that the chemical potential is 
H = 0. For numerical purposes and in all plots we set the fundamental monomer length as 
a = 1, and choose in Eq. (m) c = u/p = 1 to give a ~ 0.59(7q. The values of = 1/Re(fcq 



and 



Xq = 1/Im(fcq) correspond to the exponential decay and oscillation lengths of the g-modes, 
respectively. For fixed Xj £.q decreases and Xq increases with increasing q. Similar behavior 
was found by Petera and Muthukumar , using a different free energy functional ||] . Close 
to the ODT (but within the range of validity of the model), and for g-modes such that q > qo 
we find finite ^q and Xq ~ (xc — x)^^^^- However, all g-modes in the band < q < qo are 
equally "active", i.e., these modes decay to zero very slowly in the vicinity of the ODT as 
y ^ oo: ~ (xc — x)~^^^ is finite. Therefore, the propagation of the surface imprint 

(pattern) of g-modes with q < qo into the bulk can persist to long distances, in contrast to 
surface patterns with q > qo which persist only close to the surface. The q = qo mode has 
both lengths ^q, A, diverging as (xc - x)~^''^ for X ^ Xc- 



The boundary conditions at y = for 5(f}q are 

50q(O) - 

aq/h + 2TMq{'^)/h+{ql-q^ 



-{ql-q^)5<j}q{0) = 



(6) 
(7) 
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The amplitude Aq is found to be Aq = —aq [At, + 2lm{kq)\/Thj . Thus, in cases where the 
surface orders at the same temperature as the bulk, ~ 0, the copolymer response diverges 
upon approaching the critical point as {Nxc ~ ^x) 

With our method, any two-dimensional chemical pattern a{x, z) can be modeled. For 
surface feature size larger than do, the characteristic copolymer length, the melt can propagate 
the chemical surface pattern into the bulk. This is clearly seen in Fig. 1, where the surface 
pattern separating A regions and B regions is chosen arbitrarily to have the shape of the 
letter 'A' of size 15do x 15do- The order propagating perpendicular to the surface. Fig. 1, 
also induces order in a parallel direction, where lamellae appear oriented along the long edges 
of the letter. This ordering decays exponentially in the lateral direction, with decay length 
determined by the surface g-modes. The shading of the contours in (b) is approximately 
opposite to that of (a) and (c) indicating a reversal of the pattern; this is expected for two 
planes separated by a half-integer number of do, a distance for which an A<->B interchange of 
monomers occurs. The distance from the surface at which the pattern completely fades out 
depends on the surface Fourier components: each q-mode decays as dictated by Eq. (|5|), and 
the distant image is a superposition of all modes. In (c) the pattern is almost washed out at 
a distance 12cio from the surface. Figure 2 (a) is similar to Fig. 1 (b) (with y = 0.5do), but 
we scaled down the 'A' pattern so that its size is 4do x Ad^. For surface feature size smaller 
than do, the melt cannot follow the surface pattern, and the morphology is blurred even very 
close to the surface. As one goes deeper into the disordered phase (lowering x), the lamellar 
features are less apparent, as is seen in Fig. 2 (b) for which — 9-5 and all other parameter 
are taken as in Fig. 1 (b). 

Our treatment can be generalized to handle a thin BCP film of thickness 2L confined 
between two chemically patterned surfaces. The two corresponding surface patterns (located 
at y = ±L) are defined by two surface fields (J^{x, z) — ct^ exp [i {qxX + Qzz)]. The form 
of the response function S(j) is required now to satisfy four boundary conditions (similar to 
Eqs. (^ and (0)), two on each surface. For very large separation 2L, the melt orders close to 
the y = ±L surfaces, while the middle of the film, y « 0, is in its disordered state. As the 
inter-surface separation decreases, the ordered layers close to the two surfaces start to overlap. 

We briefly mention results for a thin-fllm BCP below its bulk ODT temperature |^^, where 
the two surface fields cr^'^ at y = ±i are taken to be uniform. Their values can be positive 
or negative depending on their preference to the A and B blocks. The free energy Eq. (||) 
expresses a correction to a lamellar phase perpendicular to the surface. The correction field 
S(f>{y) has the form Re[Ao exp(— fcoj/) -I- i?o exp(— fcg?/) + const.]. A similar perturbed lamellar 
phase was found in the strong or intermediate segregation regimes However, close to 

the ODT the lamellar deformation near the surface changes the free energy considerably. This 
is seen in Fig. 3 (a) for the symmetric case, cr+ = . The B-monomers (in black) are attracted 
to both surfaces. The oscillatory excess free energy AF < is shown in (b) as a function 
of the separation 2L between the surfaces, indicating that the deformation d(f) indeed lowers 
the total free energy. In (c) a semi-infinite system is shown, with one sinusoidally patterned 
surface of period d at ?/ = 0. The lamellae appear tilted with an angle 6 = arcsin(do/d) with 
respect to the surface, in agreement with Ref . . Close to the surface the deformation may 
be very large. 

In summary, we show a simple method to find the surface effects in confined diblock 
copolymers. The free energy is expanded to second order around a bulk (ordered or disordered) 
phase. Above the ODT we obtain a master equation for the response of copolymer melt to 
an arbitrary two-dimensional surface pattern. These two-dimensional patterns have not been 
considered previously ||l^ and show a rich behavior as the BCP order parameter replicates (to 
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some extent) the surface pattern. In order to transfer a pattern from the surface to the bulk, 
or from one surface to another via a BCP film, it is important to control the thermodynamic 
conditions, and in particular, how close we are to the ODT. Two opposite trends should be 
taken into account. On one hand, far from the ODT, the required surface fields may be 
large. On the other hand, close to the ODT, high g-modes decay much faster than low q- 
modes, resulting in a distorted pattern. Below the ODT we find the profile and energy of 
the distorted perpendicular lamellar phase as function of the surface fields and film thickness. 
When the surfaces are chemically modulated tilted lamellar phases are stabilized . 

Surface behavior of BCP or thin films of BCP confined between two surfaces may be 
used in many applications. For example, surface patterns can be used to generate a desired 
morphology inside a BCP film. Then, one of the two components is removed leaving behind 
a pre-designed porous material with controlled morphology. 
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Fig. 1 ~ Contour plots of the diblock copolymer order parameter depicting slices parallel to the y = 
surface, (a) The surface pattern at y = 0, chosen in the shape of the 'A' letter of size 15do x 15do. 
Inside the letter 'A' a — 1, while outside cr = 0. In (b) the BCP morphology is shown for y = 0.5do, 
and in (c) for y = 12do. For planes separated by a half integer number of do, an A<-^B interchange of 
monomers occurs (compare (b) to (a) and (c)). The order propagating in the perpendicular direction 
induces order in a direction parallel to the surface. This is clearly seen in (b) where lamellae form 
parallel to the edges of the 'A' letter. The Flory parameter is ~ 10-3 and seven gray scales 
are used to maximize the image contrast. All lengths are scaled by the lamellar period do, and the 
parameters chosen in Eq. (|^) are: A'^ = 1000, = ^Na^, a = c = u/p — 1 and = 0.03. These 
values are used in Figs. 2-3 as well. 




Fig. 2 - Contour plots as in Fig. 1 (b) {y = 0.5do) 
of 4do X 4do. Note that the surface cannot induce 
or smaller. In (b) the temperature is higher, xN = 
less prominent than in Fig. 1 (b). 
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, but in (a) the surface pattern is reduced to a size 
a bulk ordering when its pattern size is close to do 
= 9.5, and the lamellar features along the letter are 
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Fig. 3 - (a) Perpendicular lamellar phase between two (top and bottom) surfaces below the ODT 
(X-'V = 10.7). The two surfaces at y = zbL have equal surface fields, cr"'" = a~ = 0.07, preferring the 
B-block (in black) . (b) The oscillatory character of the correction free energy AF < (Eq. ^ of the 
symmetric system in (a), as a function of surface separation 2L. AF is divided by the bulk lamellar 
free energy Ft < 0, and it decays to zero for L — > oo. (c) Tilted lamellar phase occurring for a single 
sinusoidally patterned surface of period d > do at y = 0. The lamellae gain interfacial energy by 
overlapping with the surface pattern. The Flory parameter is =11 and the surface periodicity is 
chosen as d/do = 2 giving a tilt angle of 6 = 30°. 



